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Abstract. We study the homogenization of a Schrodinger equation in a locally 
periodic medium. For the time and space scaling of semi-classical analysis we 
consider well-prepared initial data that are concentrated near a stationary point 
(with respect to both space and phase) of the energy, i.e. the Bloch cell eigen¬ 
value. We show that there exists a localized solution which is asymptotically 
given as the product of a Bloch wave and of the solution of an homogenized 
Schrodinger equation with quadratic potential. 
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1. INTRODUCTION 

We study the homogenization of the following Schrodinger equation 

( 1 , 1 ) 1 ("■;) ° ^ 

[ Us{0,x) = u'^{x) in 

where the unknown Us{t,x) is a complex-valued function. The coefficients A{x,y) and 
c{x, y) are real and sufficiently smooth bounded functions defined for x G (the macro¬ 
scopic variable) and y G (the microscopic variable in the unit torus). The period £ 
is a small positive parameter which is intended to go to zero. Furthermore the matrix 
A is symmetric, uniformly positive definite. Of course the usual Schrodinger equation is 
recovered when A = Id but, since there is no additional difficulty, we keep the general 
form of equation fll.ll) in the sequel (which can be interpreted as introducing a non flat 
locally periodic metric). 

The scaling of fll.ip is that of semi-classical analysis (see e.g. [S], [S], [ID], [H], [T^ - 

ra. la. [18] . [T9]h if the period is rescaled to 1, it amounts to look at large, time and 
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space, variables of order At least in the case when A = Id and c(x, y) = Co(x) + Ci(?/), 
there is a well-known theory for the asymptotic limit of fll.ip when e goes to zero. By 
using WKB asymptotic expansion or the notion of semi-classical measures (or Wigner 
transforms) the homogenized problem is in some sense the Liouville transport equation 
for a classical particle which is the limit of the wave function u^. In other words, for an 
initial data living in the n-th Bloch band and under some technical assumptions on the 
Bloch spectral cell problem fll.4l) . the semi-classical limit of fll.ip is given by the dynamic 
of the following Hamiltonian system in the phase space (x, 6) G x 


( x = VeXnix,e) 

\ 0 = -V,A„(x,0) 

where the Hamiltonian Xn{x,6) is precisely the n-th Bloch eigenvalue of fll.dp (see [H], 
pH] . [TT] . [12], [13], [H], [18], [19] for more details). 

Our approach to fll.ip is different since we consider special initial data that are monochro¬ 
matic, have zero group velocity and zero applied force. Namely the initial data is concen¬ 
trating at a point (x'^, 9"’) of the phase space where VeA„(x”, 6^) = VxXn{x"', 0”) = 0. In 
such a case, the previous Hamiltonian system fll.2|) degenerates (its solution is constant) 
and is unable to describe the precise dynamic of the wave function u^. We exhibit another 
limit problem which is again a Schrodinger equation with quadratic potential. In other 
words we build a sequence of approximate solutions of fll.ip which are the product of a 
Bloch wave and of the solution of an homogenized Schrodinger equation. Furthermore, if 
the full Hessian tensor of the Bloch eigenvalue An(x, 9) is positive dehnite at (x”, 6*”), we 
prove that all the eigenfunctions of an homogenized Schrodinger equation are exponen¬ 
tially decreasing at infinity. In other words, we exhibit a localization phenomenon for fll.ip 
since we build a sequence of approximate solutions that decay exponentially fast away 
from x"'. The root of this localization phenomenon is the macroscopic modulation (i.e. 
with respect to x) of the periodic coefficients which is similar in spirit to the randomness 
that causes Anderson’s localization (see [9] and references therein). 

Let us describe more precisely the type of well-prepared initial data that we consider. 
For a given point (x*^, 0”') G x and a given function G we take 

(1.3) «;(!) =V>„(x", 

where 'ipn{x, y, 9) is a so-called Bloch eigenfunction, solution of the following Bloch spectral 
cell equation 


(1.4) 


(divy -f 2i7r6')(A(x, i/)(Vy 2i7r6')i/>„) -1- c(x, y) = A„(x, 9)i>n 


in T^, 
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corresponding to the n-th eigenvalue or energy level A„. The Bloch wave V'n is periodic 

x-x" 

concentrated around with a support of asymptotic size y/e. The Bloch frequency 
Qn g 'jN^ -j-pg localization point x” G and the energy level n are chosen such that 
Xn{x^,0^) is simple and VxAn(x”, 6*"') = V6)A„(x”', 6*"') = 0. 

Our main result (Theorem 13.2p shows that the solution of fll.ip is approximately given 
by 


^ means that the initial data is 


with respect to y but is not periodic, so 


(1.5) 


Ue{t,x) K. 


9")f g, 9".^ / X - X- 

* ^ ^ Vit, 


V~e 


where v is the unique solution of the homogenized Schrodinger equation 


( 1 . 6 ) 


i— -div (A* Vn) + div{vB*z) + c*v + vD*z ■ z 

at 

^(0, z) = v^{z) 


0 in X M+ 
in 


where c* is a constant coefficient and A*,B*,D* are constant matrices dehned by 
Al* = ^VeV,A4x^r), B* = ^V,V,A„(x^0"), D* = ^V,V,A„(x^n • 

OTT^ /ZTT / 


In Proposition 13.41 we show that the homogenized problem fll.bp is well-posed since the 
underlying operator is self-adjoint. Furthermore, under the additional assumption that 
the Hessian tensor VVA„(x'^, 0”') (with respect to both variables x and 6) is positive 
dehnite, we prove that fll.bp admits a countable number of eigenvalues and eigenfunctions 
which all decay exponentially at inhnity (see Proposition 13.5p . In such a case, formula 
fll.bp dehnes a family of approximate (exponentially) localized solutions of fll.ip . 

Let us indicate that the case of the hrst eigenvalue (ground state) n = 1 with 6^ = 0 
was already studied in [3] (for the spectral problem rather than the evolution equation). 
The case of purely periodic coefficients (i.e. that depend only on y and not on x) is com¬ 
pletely different and was studied in |1]. Indeed, in this latter case there is no localization 
effect and one proves that, for a longer time scale (of order with respect to fll.lD L 
the homogenized limit is again a Schrodinger equation without the drift and quadratic 
potential in fll.bp . 


2. Preliminaries 

In the present section we give our main assumptions, set some notation and a few pre¬ 
liminary results needed in the proof of the main results in Section [3l 

We hrst assume that the coefficients Aij{x,y) and c{x,y) are real, bounded, and 
Caratheodory functions (measurable with respect to y and continuous in x), which are 
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periodic with respect to y. In other words, they belong to Cb (M^; . Further¬ 

more, the tensor A{x, y) is symmetric uniformly coercive. Under these assumptions, it is 
well-known that, for any values of the parameters 6 G and x G M^, the cell problem 
fll.411 dehnes a compact self-adjoint operator on which admits a countable se¬ 

quence of real increasing eigenvalues {A„(a;, 9)}n>i (repeated with their multiplicity) with 
corresponding eigenfunctions {'ijjn{x,9,y)}n>i normalized by 

\\i^n{x,9,-)\\L2(^jN) = 1 . 

Our main assumptions are: 

Hypothesis HI. There exist x"' G and 9^ G such that 


( 2 . 1 ) 


{i) Xn{x^, 9^) is a simple eigenvalue, 

(n) {x'^,9^) is a critical point of \n{x,9),i.e. VxXn{x"',9"') = VeA„(a;"^, 6*”) = 0. 


Hypothesis H2. The coefficients A{x,y) and c{x,y) are of class with respect to the 
variable a: in a neighborhood of a: = x". 


Then we set: 


(9/1 (9^/1 

^hh{y) ■= ^:-{.x"',y), A,ih{y) •= {x^,y), for i,h = i, 


.N. 


v ' ^ J ' ^.Llt O 

OXh OXiOXh 

Similar notation is used to denote the derivatives of the function c with respect to the 
x-variable. With an abuse of notation we further set 


A{y) ■.= A{x\y), A„:=A„(x",n, My) ■= i^n{x\y,9^), 

and analogous notation holds for all derivatives of xfjn and A„ with respect to the x-variable 
and the 0-variable evaluated at x = x"" and 9 = 9"'. Without loss of generality we will 
assume in the sequel that x" = 0. 

Notation. For any function p{y) dehned on we set 

p%z) := p{z/M) 

where z := ^/ey = xjs/e. In the sequel the symbols divj, and Vy will stand for the 
divergence and gradient operators which act with respect to the ?/-variable while div 
and V will indicate the divergence and gradient operators which act with respect to the 
2 ;-variable. Finally throughout this paper the Einstein summation convention is used. 

Under assumption fl2.1ll -fil it is a classical matter to prove that the n-th eigencouple of 
fll.41) is smooth with respect to the variable 0 in a neighborhood oi 9 = 9" (see [IB]) and 
has the same differentiability property as the coefficients with respect to the variable x. 
Introducing the unbounded operator A„(x,6*) dehned on L'^(T^) by 

A„(x, 9)'ip = —(divy -|- 2i7i9) (^A{x, y)(Vy + 2i7i9)'ip^ + c(x, y)'ip — A„(x, 9)'ip, 
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it is easy to differentiate fll.dp . Denoting by (efc)i<A:<Ar the canonical basis of , the hrst 
derivatives satisfy 


( 2 . 2 ) 


An{x,9) 


djjn 

d6k 


2inekA{x, |/)(Vj/ + 2i7T6)i>n 

dX 

+ (divj^ + 2me) {A{x, y)2mek'ipn) + 9)'4)n , 

CWk 


(2.3) 


Ayi (x, 9 ) 


djJn 

dx. 


/ d A 

(divy + 21719) ^)C^y + 2,i7i9)'il)y 

f)r r)X 


Similar formulas hold for second order derivatives. By integrating the cell equations for 
the second order derivatives against ipn we obtain the following formulas that will be 
useful in the sequel (their proofs are safely left to the reader). 

Lemma 2.1. Assume that assumptions HI and H2 hold true. Then the following equal¬ 
ities hold: 


(2.4) 




^ ^A,^k{Vy + 2z7rd-)^ • (V, - 2^7r0-)^„ + ci,, 

OUk Ot>k 


27ri . 


+ 


IjN L 


^^|Jn 


Al^hGki’n ■ (Vj^ - 2i7r6'”)V’n + ■ (Vj/ - 2i7r6'”)'0n 

dxh 


eki’nAi^h ■ (Vj^ + 2i7T9'^)i!n + eki>nA ■ {Vy + 2i'K9'^) 


JjN L 

1 d^X^ 
2i'K dxhd9k 


dy 

dy 

dy 


= 0 , 


(2.5) 


JjN L 

+ 


02 \ . 

A2,lh{'^y + 2i'K9'^)'lpn ' (Vj^ - 2i7r6'’")'0„ + ( C2,lh - IV'nP dy 


IjN L 


+ 


IjN L 


24i,;,(V, + 2*7rd")^ • (V, - 2md")^„ + 
Ai/V, + 2^vr0")^ • (V, - 2i-K9^)fij^ + Ci,, 


dy 

dy = 0, 
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( 2 . 6 ) 


?TiV 


2i7iekA{y)(yy + - {^{y)2mek^^^ ) (Vy - 2i7r6'”)V>„ dy 


d'lpn 


+ 


IjN 


/JN 

d^X 


2i7ieiA{y){Vy + 2z7r6''^)^^V’n - ( A{y)2mei^^ ) (Vy - 2i7r6'”)V'„ dy 


dA. 


dOk 


ATT‘^ekA{y)ei\'4)n\^ + ATT‘^eiA{y)ek\'4)n\^ dy 


dOk 


^ {nn 


(0”) = 0. 


dOidOk 

We now give the variational formulations of the above cell problems, rescaled at size £. 

Lemma 2.2. Assume that assumptions HI and H2 hold true and let (p(^) he a smooth 
compactly supported function defined from into C. Then the following equalities hold: 


(2.7) 


7l^(Vy + 2z7r0"X • (^/iV - 2md'^)^{z) + {d - X^fj^T 


dz = 0, 


( 2 . 8 ) 


V(V„ + 2»>r9")^ . (ViV - 2^„0'')g, + (c' - K)^J 


dz 


2'Kiek ■ 24^(Vy + 2f7r6*”)'0^(p + A^ 2'Kiek'f’n ' (\/^^ ~ 2i7r6*"')(p dz = 0 , 


(2.9) 


A’(Vy + 2 !ir 9 ")hS . (^v - 2!ir9”)<^ + (c' - A')^<^ dz 


dxh 


dxh 


+ 


• (\/iV - 2z7rd")v7 + 


dz = 0. 


Proof. Formula fl2.7p follows straightforwardly from equation fll.4p while fl2.8l) - fl2.9p are 
consequences of fl2.2p - fl2.3l) . □ 

Finally we recall the notion of two-scale convergence introduced in [T], [T7] (that will 
be used with 6 = y/e). 

Proposition 2.3. Let fs be a sequence uniformly bounded in L‘^{W^). 

(1) There exists a subsequence, still denoted by fs, and a limit fo{x, y) G x T-^) 

such that fs two-scale converges weakly to fo in the sense that 

lim / fs{x)(j){x,x/6)dx= / fo{x,y)(j){x,y) dx dy 

< 5^0 J^N J^N JjN 

for all functions 0(x, y) G C'(T^)). 
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(2) Assume further that fs two-seale converges weakly to fo and that 

||/5||l2(R^) = ||/o||l2(R^xT^)- 

Then fs is said to two-scale converge strongly to its limit fo in the sense that, if 
fo is smooth enough, e.g. fo G Lf (M'^; C'(T'^)), we have 

lim [ \fs{x) - fo{x,x/6)f dx = 0. 

J^N 

(3) Assume that SV fs is also uniformly bounded in . Then there exists a 

subsequence, still denoted by fs, and a limit fo{x,y) G L^(R^; H^(T^)) such that 
fs two-scale converges to fo{,x,y) and SVfs two-scale converges to Vyfo{x,y). 


3. Main results 


We begin by recalling the usual a priori estimates for the solution of the Schrodinger 
equation dUD which hold true since the coefficients are real. They are obtained by 
multiplying the equation successively by and and integrating by parts. 


Lemma 3.1. There exists C > 0 independent of e such that the solution of fll.ip satisfies 


I |W| |l°°(R+;L2(rJV)) — I |m°| |2,2(RiV) , 

^11 Vrtel |Loo(R+.i2(RiV)) < C'f I |L2(RiV) + e:| | Vm°| |i2(RjV) 


Theorem 3.2. Assume that assumptions HI and H2 hold true and that the initial data 
is of the form fll.31) . Then the solution of fll.ip can be written as 


(3.1) 


Ue{t,x) = e^^e^^-—vAt 


X — X 




where Vs(t,z) two-scale converges strongly to 'ifn{y)v(t, z). 


(3.2) 



Ve{t,z) - fjn 



i.e. 

2 

dz 


0 , 


uniformly on compact time intervals andv is the unique solution of the homogenized 

Schrodinger equation 


(3,3) 


3v 

i— -div (y4* Vn) + (\.yv{vB*z) + c*v + vD*z ■ z = in x R+ 

at 

t( 0 , z) = v^{z) in 


dV 


where 


= -^^eVe\n{x^,eA , B* = ^V,V,A„(a;^0"), D* = 0”), 

StT^ IlTT 2 
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and c* is given by 


c = 


I'lN L 


diJj d%b — 

A{Vy+2i'K9'^)'4)n-^^ek-A{Vy-2i'K9'^)—^-'4)n ek-Ai^k{'^y-2i7i9^)ipn-'ilJnek dy . 

oxk dxk 


Remark 3.3. Notice that even if the tensor A* might be non-coercive, the homogenized 
problem fl3.3p is well posed. Indeed the operator A* : —)■ L^(R^) dehned by 

(3.4) A*Lf = —div (A* V</9) + div{ipB*z) + c*cp + ipD*z ■ z 

is self-adjoint (see Proposition 13.4p and therefore by using semi-group theory (see e.g. [B] 
or Chapter X in EDI). one can show that there exists a unique solution in C(R’''; L^(R'^)), 
although it may not belong to L^(R+; if^(R^)). 

The next result establishes the conservation of the L^-norm for the solution v of the 
homogenized equation fl3.3|) and the self-adjointness of the operator A*. 

Proposition 3.4. Let V e C(R’''; L^(R^)) be solution to fl3.3p . Then 

(3.5) ||'i^(C OIIl^cr^) = ll'^°llL2(RiV) VtGR’*'. 

Moreover the operator A* defined in 03.41) is self-adjoint. 

Proof. We multiply the equation 03.3p by v and take the imaginary part to obtain 

(3.6) f \v\‘^dz = lm( f vB*z ■ Vv — G 

2 dt Jrjv \Jrn 

After integrating by parts one hnds that the right hand side of 03.6p equals 


—tr R*-|-Imc*) / Inpdz 
2i / J^N 

and therefore 03.5p is proved as soon as we show that 


(3.7) 


—tr B* + Imc* = 0 . 
2i 


In order to do this we hrst rewrite the coefficients c* and B* in a suitable form. Denoting 
by (•, •) the Hermitian inner product in L^(T^) and using equation 02.2p we write 


(3.8) 

while by equations 02.2p 


c* = - [ Ai^k{yy-2in9'^)fin-'finekdy, 

2x71 d9k oxk 


it follows that 


(3,9) 


1 ^_p pp _p, pp pp, 

2i7r dxhd9k 2m ” d9k ’ dxh 2m ” dxh ’ d9k 


22 lm / Ap/,(Vy - 2m9'^)fin ■ dy ■ 
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By formulae fl3.8l) - fl3.9p it is readily seen that equality fl3.7p holds true. 

In order to prove the self-adjointness of the operator A*, one hrst checks that A* is 
symmetric, which easily follows by fl3.7p and the fact that B* = —5*, and then observes 
that up to addition of a multiple of the identity the operator A* is monotone (see e.g. [7], 
Chapter Vll). □ 


In the next proposition we will denote by VVAn the Hessian matrix of the function 
Xn{x,9) evaluated at the point {x'^,9'^), namely 


VVAfi 


^xXn ^8^xXn \ / n 

VeV.A, VgVeXn J > 


Proposition 3.5. Assume that the matrix VVA„ is positive definite. Then there exists 
an orthonormal basis {(pn}n>i of eigenfunctions of A*; moreover for eaehn there exists a 
real constant 7 n > 0 such that 

(3.10) G L\R^) . 


Proof. Up to shifting the spectrum of the operator A*, we may assume that Re(c*) = 0. 
In order to prove the existence of an orthonormal basis of eigenfnnctions we introduce the 
inverse operator of A*, denoted by G* 

G* : L^iR^) ^ L^iR^) 

/ ^ (p nniqne solntion in H^iR^) of 

(3.11) AV = / in 

and we show that G* is compact. Indeed mnltiplication of fl3.1ip by yields 

(3.12) / [A*V‘p-Vf) — iB*\vsi{ipz-Vifi) + D*z ■ z\ip\^]dz = f ffidz. 

Upon dehning the 2iV-dimensional vector-valned fnnction $ 

■ ) 

we rewrite fl3.12p in agreement with this block notation 

^ VVAnd) ■'^dz = 


Stt^ 


fip dz. 


By the positivity assnmption on the matrix VVA^ it follows that there exists a positive 
constant Cq snch that 


Co(^||Vv?||^2(]RjV) + |k<7’|lL2(]RjV) j < ||/||L2(RiV)||(y9||i2(RiV) , 

which implies by a standard argument 

I IT^I lL2(KiV) + I I V(p| |j;^2(RiV) + I lL2(RiV) < ^ \\ f \\ j^2 fT^N) , 
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from which we deduce the compactness of G* in L^(M'^)-strong. Thus there exists an 
infinite countable number of eigenvalues for A*. 

We are left to prove the exponential decay of the eigenfunctions (this fact is quite 
standard, see e.g. 0 )- Let ipn be an eigenfunction and let be the associated eigenvalue 

(3.13) A*ipn = aniPn- 


Let > 0 and p G C'°°(M) be a real function such that 0 < p < 1, p(s) = 0 for s < i?o 
and p{s) = 1 for s > i?o + 1 and for every positive integer k define pk G C'°°(M^) in the 
following way 

Pk{z) := p{\z\ - k). 

We now multiply 03.131) by ’^npi to got 


pI {A*Vp>n ■ Vpr. 


■ Vp>n) + D*Z ■ z\p>n\^ - an\p>n?) dz = 


(3.14) 


{pk\p)n?B*Z ■ Vpk - 2pk iPnA*Vlfn ' Vpfc) dz . 


Next remark that since the left hand side of 03.141) is real the right hand side must be 
also real and therefore it is equal to 


(3.15) 


-2pk Re(p„A* V<p„) • Vpk dz. 


Let Bk denote the ball of radius Rq + k and center z = 0 and observe that the support 
of Vpk is contained in Bk+i \ Bk. Then putting up together 03.14p and 03.15p and using 
again the positive definiteness of the matrix VVA„ we obtain for Rq sufficiently large 
(\/^ > an does the job) 


I (^1 I^Pnl I 

where ci is a positive constant independent of k. Thus we deduce that 


(3.16) 


Ha 




Cl 


. 1 + Cl 

Upon defining a positive constant 70 > 0 by 


11 A' 


n\\m{RN\Bo) 


Cl 


1 + Cl 


_ g-27o(fc+Ro) 


it is finally seen that 03.16P implies the estimate 03.1Up for any exponent 0 < 7 ^ < 70 . □ 

Proof of Theorem 13.21 We rescale the space variable by introducing 

X 


z = 




and define the sequence v^. by 

(3.17) 


/, \ —jhnl 2?T /, \ 

Ve{t,z) := e s e ^ ® Ue{t,x) . 
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By the a priori estimates of Lemma [3.II it follows that v^it, z) satishes 

ll'f^£||L°°(R+;L2(RiV)) + \/£||Vi;£||j;,oo(R+.i2(RjV)) < C, 

and applying the compactness of two-scale convergence (see Proposition 12.31) . np to a 
snbsequence, there exists a limit v*(t,z,y) G (M’*' x H^(T^)) snch that and 
y/eVve two-scale converge to v* and 'VyV*, respectively. Similarly, by dehnition of the 
initial data, v^{0,z) two-scale converges to ' 0 „(i/)n°( 2 ;). 

Althongh Vg is the unknown which will pass to the limit in the sequel, it is simpler to 
write an equation for another function, namely 

(3.18) We{t, z) := ng(t, z) = e~"~^ Mg(f, x). 


By fl3.18p it follows that 

(3.19) 


VtCg = e -|- 2i7i-^jVs , 


and it can be checked that the new unknown solves the following equation 

(3.20) 

dw, 


div[24 (^\/ez, z/\/e) VtCg] H— \c{\fez^zl\fe) — Xn]ws = 0 in x M+ 

CJTi S 


m 


f,N 


Wg(0,2;) = Mg(v^^) 

where the differential operators div and V act with respect to the new variable z. 
First step. We multiply the equation fl3.20p by the complex conjugate of 


, / Z \ 2lTT^ 

£(p(t,z,-^y ^ 


where 0(s, z, y) is a smooth test function dehned on M’*' x R^ x T^, with compact support 
in M"*" X R^. Since this test function has compact support (hxed with respect to e), the 
effect of the non-periodic variable in the coefficients is negligible for sufficiently small 
e. Therefore we can replace the value of each coefficient at {\/ez, z/y/e) by its Taylor 
expansion of order two about the point (0,z/i/e). Integrating by parts and using fl3.18p 
and fl3.19p yields 


p + OO ( 


IS 




dtdz — ie / Ve{0, z)4>(o, z,dz 

Jrn V ^ye/ 


Jo Jm.^ 

P+OOP 

+ [A^ + Al i^ '/ezh + 1 ^ 2 eziZh + 0 (e)] (i/eV -h 2m6'^)v^- {^/e'V - 2i7r6^)(j)‘^ dz dt 

Jo Jm.^ 

P+OOP 

+ / [c" + cl^hV^Zh + \cl^ih ^ziZh + o(e) - A„] Ve4>^ dz dt = 0. 

Jo Jr’^ 

Passing to the two-scale limit we get the variational formulation of 

— (divj^ -|- 2^6*") (^A{y){'Vy + 2i7r6*”)n* j -|- c{y)v* = in T'^. 
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The simplicity of implies that there exists a scalar function v{t,z) G (M"*" x 
such that 


(3.21) 


v*{t,z,y) = v{t,z)^pn{y)■ 


Second step. We multiply fl3.20p by the complex conjugate of 




^ f 1 d(j) ^ ^ 

k=l 


where 0(t, z) is a smooth test function with compact support in M’*' x We first look 
at those terms of the equation involving time derivatives: 


^H-OO / 


dw^ 


(3.22) / / p dt dz = 


IQ Jr^ 

p+oop 

'0 Jr^ 


dt 


-IVp 


N 




dt 


k=l 


1 d'lpn d'^4> d'lpn d4> 
T Zk 


2i'K dOk dtdzk dxk dt 


dt dz 


i / VeiQ.z) 
Jr^ 


N 

C0(O,^) + 

k=l 


1 dip^ d(j) , s dil)i . 

2171 dOk dzk dxk 


dz. 


Recalling the normalization Jjpjv dy = 1, we find that the two-scale limit of the term 
on the left hand side of fl3.22p is given by the expression 


(3.23) 


(* + 00 / 


dcf) 


— i / v-^dzdt — i / v^(j)(0, z) dz. 

Jo Jm^ dt J^N 


We further decompose as follows 


= + z with 



N 


k=l 


d'iPn 

dxk 


((>{t,z)ek. 


Getting rid of all terms multiplied by o{e) and taking into account fl3.18p and fl3.19p we 
next pass to the limit in the remaining terms of fl3.2Up multiplied by The computation 
is similar to |1] but it involves new terms since ipn and its derivatives also depend on x. 
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We first look at those terms which are of zero order with respect to z, namely 


(3.24) [ \a^Vw, • + ^2) + -(c" - Xn)We^l 


>0 JR^ 

t' + OOp 


dz dt 


-AA^/eV + 2i7r9'^)vs ■ (V^ - 2i7r6'")'0^(?!) + -(c^ - XnWnVecf) 
0 L£ V / £ 


dz dt 


2i7r 


^+oo / 


'0 Jr^ 


A^(^^/eV + 2i7re^y, ■ {Vy - 2^6^) 


We 

, 1 ^ ^ d^j^dct) 

^/e 89 k dzki 


dipn d(j) 
d9k dzk 


dz dt 


+ f f + 2z7r6'"'')ne • 

Jo Jr^ -s/s \ / 

+ f f -7^4l^(x/iV + 227r0’^)n,- 

Jo Jr^ 27rz V / 


'0 Jr^ 


d9k dzk 

djj^ 

Vp ■ ^r^<pek dzdt. 


dxi 


d(j) 

Using equation fl2.7p with ip = v^cp and equation fl2.8p with ip = v^— — we rewrite the hrst 

8zk 

two integrals in the right hand side of fl3.24p as follows 


f [ — ^A'^{'Vy — 2i7i9'^)il)l^ ■ VeVJfdzdt 

0 Jr^ V ^ 

+ n \W(^, - 


'0 JR^ 


.2171 


1 


- -j=A^il)l^ek ■ (-v/eV + 2m9'' 

ye 


d(j) \ 

' 8 Zk) 


dz dt. 


Combining the above terms with the other terms in (13.241) and passing to the two-scale 
limit in fl3.24p yields 


(3.25) 


n / 

'0 2rJV JtJv 

^ /7 J 

Jo Jr^ Jt^ 


—A'lpnWy - 2i7r6'”)^2i _ + 2i7r6'”)'0n - ^dlV'npefc 

L2z7r 89 k 2nr 89 k 

dcj) 

vV —— dy dz dt 
8zk 


A{Vy + 2m9^)i)r^ 


8J)r. 

8xk 


■v(j)ek dy dz dt. 
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By equation fl2.6p it can be seen that the hrst integral of 03.251) equals 

/*+oo /* 


(3.26) 


A*'Vv'V(j)dz dt. 


Jo 

We now focus on those terms which are linear in z: 


^H-OO / 


'0 Jr^ 


■ (V'h'z) + -{d - \n)w,^z + Al^^^ezuVw, ■ + 'h') 




cUZkW^^l 


dz dt 


f+cor 

Jo Jr^ 

p+oop 

'0 Jr^ 

P+OO P 


1 / \ - 1 - 
A^(^/iV + 2i7r6'”)n£- {Vy - (pZk + —{c^ - Xn)vs^r:^ 4>Zk 


We 




dxi 


dz dt 


We 


-7^W,kiVe^ + 2i7ie^)ve-Wy - 2iTie^Wn Wk + -^WWn^Zk 


d^j', 


0 Jr^ 

'•H-oo 


A%W^ + 2i7ieWe ■ + AlkiVe^ + 2me^)ve • 


dz dt 
dz dt 


2in J 

/* + OD 


dx_ 

+ 2i7ie'^)Ve ■ Wy - 2*716'’")^^ Zh + cl hVe^^^ Zh 


'0 Jr^ 


0 JR^ 

y/e^i hiVe'^ + 2i7i9We 


dOk dzk ’ d6k dzk 

1 dW VT W dW 7 \ j 

-V- -h — (p^kjZh dzdt. 


dz dt 


2iTi dOk dzk dxk 


(3.27) 


By equation 02.9p with ip = v^cjyzk it can be seen that the sum of the hrst two integrals in 
the right hand side of 03.271) gives 

(3.28) - r f A%Vy-2t7re^)^-v,V{Wk) + AlkWy-2t7reWn-Ve^{Wk))dzdt. 
Jo Jr^ eJXk ’ / 

Therefore passing to the two-scale limit in 03.27P we hnd 


(3.29) 


n / 

Jo Jr^ Jt^ 

n I . 

p+oo 

Jo Jr^ JfN I 

± n 

2m 


A{S/y — 2m9 


n^Wn 

dxk 
djjr 
dxk 


vi)n(t) Cfc Ai^kWy - 2m9Wn ■ v^JnCf) Ck dy dz dt 

dy dz dt 


A{Vy - 2m9'^)—J- ■ vipnZkVcf) + Ai^kWy - 2m9^)i!n ■ vipnZkVcf) 


d%b — _ _i 

AiVy + 2m9'^)'ilJn ■ v—-22zkV(f) + Ai^kWy + 2m9"')'iljn ■ vJjnZkVcj) dy dz dt 

OX}- 


h JR^ JT^ 


A.^kWy + 2*7r0-)V>„ ■ (V, - 2t7T9^)^vZkp- 

oUk ozk 


dipn d(f) 

+ Ci,hW^VZh^ 
o^k OZk J 


dy dz dt. 
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By equation 02.41) it follows that the last integral in 03.291) is equal to 


(3.30) 


n I 

P + OO p , 


L... dcf) 


Ai,hi>nek ■ (Vy - 2me'^)i)n + Mndk ■ (Vy - 2i7r6'”)-— vzh^ dydzdt 

OXh J OZk 


fo 

p-\-oo p p 


— Oib 

Al^hilnGk ■ (Vy + 2i'K9^)'llJn + ' (Vy + 2i7r6'")'0^ 


^0 , , ,, 
vzhT ,— dydzdt 
OZk 


I d K I , ,2 ^ ^ 

' / / 77 —^— ^\Vn\ vzh^ dydzdt. 

0 Jm Jt^ 2z7r dXhOOk dzk 


Next notice that the hrst and the second line of O3.30p cancel out with the second and 
the third line of 03.29p respectively and therefore 03.29p reduces to 


(3.31) 


“+ 00 / 


'0 Jr^ Jt^ 


dib 

A{Vy - 2i'K6'^)-^ ■ V'ijjn4>ek + ^i,fc(Vy - 2z7r6'”)'0n • v'lpnfpek 


dy dz dt 


/•+7 1 d^X^ ^ y y, 

Iq J^n 2m dxhddk^ dzk^'^ ^ 


Finally we consider all quadratic in z terms: 


^+OD / 


2 Jo Jrn 

+ 07 

Jo Jr^ 

- \ n 

^ Jo Jr^ 

- \ n 

^ Jo Jr^ 


'•+00 / 


Ah^eziZhVwe ■ (Vd'^ + + ch^ZiZhWe^l 


Af f.^ZkVw^ ■ (zV'If^) + k^kWeZ ■ 

ve ’ 


dz dt 
dz dt 


Z^2,lh 


'ey + 2me'^]v, 


^^(Vy-2mn0n0 + 0nV0' 
\/e 


dz dt 


^2,ih '/£ZiZh(yeV + 2me'^^Vp 


1 d(j) ^ 1 di2f,^d(j) , djjf^ 

o-^ V ^ Vt^ -h Cfc- 

2m dOk dzk 


/ / A\j^Zh[ y/eV + 2me' 

'0 ^ 

“+0O/- 




2m dOk ' dzk ' """ dxk 

ZkiWy - 2m0")^0 + ^ezk^W^ 
1 (?7n 00 


dz dt 


>0 9r^ 

/‘+°°/‘ djj^ - 

/ / c\f^ZhVeZk^^(l)dzdt 

h Jr^ 
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which give on passing to the two-scale limit 
(3.32) 

-CXD p p 

+ 2z7r6'")'0n ■ (Vy - 2z7r6'”)'0„ C2,lhiJn^^ 


1 


0 JT^ 

/*+oo p 


Al^h{^y + 2iTie'^)ipn ■ {^y - 2i7r6'”)^^ Ci^h'lpn^ J 


v(j) ziZh dy dz dt 

v(j) ZhZk dy dz dt. 


Jo Jr^ Jt^ 

Now using equation fl2.5l] we find that 03.321) reduces itself to 


P + CXD / 


(3.33) 


1 


I^N 2 dxidxh 


v(j)ZiZh dzdt. 


Summing up together 03.231) . 03.251) . 03.261) . 03.311) and 03.33P yields the weak formulation 
of 03.3p . By uniqueness of the solution of the homogenized problem 03. 3 p , we deduce that 
the entire sequence two-scale converges weakly to 'ilJn{.y)'vit,x). 

It remains to prove the strong two-scale convergence of v^. By Lemma [3.11 we have 


|'i^e(^)||L2(R^) — ||'We(^)||L2(R^) — I kel |L2(RAr) ^ | |l2(r]VxTJV) — | |n°| |i2(R]V) 


by the normalization condition of 'ijjn- From the conservation of energy of the homogenized 
equation 03.31) we have 

II'^(^)IIl2(RJV) = I |n°| |2,2(RiV), 

and thus we deduce the strong convergence from Proposition 12.31 □ 


Remark 3.6. As usual in periodic homogenization [0,0, the choice of the test function 
'kg, in the proof of Theorem 13.21 is dictated by the formal two-scale asymptotic expansion 
that can be obtained for the solution of 03.201) . namely 


We{t, z) ZZ e 




where v is the homogenized solution of 03.3p . Actually the homogenized equation that 
one gets by the asymptotic expansion method is 

dv 


(3.34) 


dt 


div (A* Vn) -|- B*'Vv ■ z + c*v + vD*z ■ z = 0 , 


which apparently differs from 03.3p by the following zero-order term 


(tr (VyV^An) - 47rlm(c*)) v . 

By virtue of 03.7p the above term vanishes, so that formulae 03.34p and 03. 3 p are equivalent. 
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